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This  contract  "Laser  and  Optical  Physics"  has  dealt  with 
problems  concerning  electromagnetic  interactions  with  matter.  On 
the  one  hand,  we  have  been  concerned  with  the  operation  of  lasers 
themselves  -  particularly  at  high  powers.  On  the  other  hand,  we 
have  been  concerned  with  laser  interactions  with  matter  as  a 
means  of  studying  the  matter  itself  and  its  polarization  carrying 
excitations. 

Laser  F ield  Distributions. 

The  design  of  lasers  requires  the  ability  to  calculate  the 
electromagnetic  field  distribution  for  a  general  configuration  of 
mirror  separation  and  radii.  For  mirrors  with  sharp  edges, 
particularly  rectangular  mirrors  such  as  those  used  in  TEA  and 
similar  lasers,  previous  calculational  methods  have  severlv  taxed 
the  facilities  of  the  worlds  largest  computers.  One  of  the 
crowning  achievements  of  the  present  contract  is  the  development 
of  a  set  of  programming  tools  that  in  combination  have  been  used 
to  calculate  electromagnetic  field  distributions  in  a  reasonable 
time  on  such  a  modest  computer  as  the  PDP-10.  These  tools 
include 
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The  use  of  horizontal  and  vertical  reflection  symmetry  to 
reduce  by  a  factor  of  four  the  core  storage  required. 

The  use  of  Fourier  transform  algorithms  that  take  advantage 
of  this  symmetry  to  reduce  computation  time  by  a  similar 
factor.  The  Fourier  transform  algorithms  avoid  certain 
instabilities  associated  with  finite  difference  algorithms.. 

The  use  of  simultaneous  integration  oi  forward  and  backward 
waves. 1  This  permits  a  reduction  of  storage  by  a  factor  M 
where  M  is  the  number  of  gain  sheets  used.  This  is 
typically  a  factor  of  6  to  10.  The  combined  saving  in 
storage  with  item  (1)  above,  is  then  a  factor  24  to  40.  The 
running  time  for  simultaneous  integration  doubles  per  loop 
between  the  mirrors,  but  the  numbers  of  loops  is  reduced  by 
more  rapid  convergence. 

A  proper  account  is  taken,  for  the  first  time,  of  the 
interference  between  forward  and  backward  waves. 2 

Propagation  from  the  sharp  and  edged  mirrors  is  handled  by  a 
continuous  Fourier  transform  (CFT)  proqram^  developed 
specifically  for  this  problem.  This  CFT  program  is  based 
upon  a  spline  fit  to  the  function  followed  by  an  analytic 
integration  of  the  B-spline  basis  functions.  With  the 
application  of  suitable  end  corrections,  we  have  achieved 
using  256  points  high  (five-figure)  accuracy  for  calculating 
Fourier  integrals  of  functions  with  jumps  or  cusps. *  The 
usual  Fast  Fourier  Transform  algorithms  with  the  same  data 
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give  much  lower  accuracy  (two  figures)  or  require  an 
enormously  high  number  of  points  to  achieve  the  same 

accuracy  as  CFT.  Thus  our  CFT  program  is  often  faster  than 
FFT,  for  the  same  accuracy,  for  functions  with 

discontinuities  or  kinks.  Moreover,  our  CFT  program  has  not 
yet  been  written  in  a  form  to  achieve  high  speed. 

It  should  be  emphasized  that  the  procedures,  described 
above,  as  most  successful  in  solving  the  nonlinear  partial 
differential  equations  of  the  electromagnetic  field  were  not 
proposed  in  the  original  contract  proposal.  Contract  proposals 
should  not  be  expected  to  provide  explicit  indications  of  the 
methods  to  b.  used  -  since  these  predictions  will  be  wrong, 
unless  one  already  knows  the  answers. 

Optical  Bistability. 

We  have  also  recognized  the  importance  of  optical 
bistability  for  optical  switches,  memories  and  amplifiers.®*®  We 
have  therefore  made  a  study  of  bistability  induced  by  nonlinear 
dispersion,  in  addition  to  the  nonlinear  absorption  usually 
concerned.  This  problem  was  noi.  even  among  those  originally 
proposed . 

Laser  Annealing. 

Another  area  of  importance  to  develop  recently  is  the  ion 
implantation  method  of  dopants  in  semiconductor  devices.  This 
precision  method  of  doping,  however,  damages  the  surface.  The 
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damage  can  be  removed  in  milliseconds,  or  less,  by  annealing  the 
surface  with  a  scanning  laser  beam.  [Laser  annealing  has  been 
the  subject  of  invited  sessions  at  the  American  Physical  Society 
and  the  Materials  Research  Society  recently.]  A  proper 
understanding  of  the  kinetics  of  annealing  requires  a  knowledge 
of  the  temperature  rise  induced  by  the  laser  beam.  Although  the 
temperature  rise  for  a  weak  laser  beam  is  a  "textbook"  problem, 
we  published  the  first  solution  to  this  problem.7  However, 
annealing  requires  a  strong  laser  beam  which  raises  the 
temperature  almost  to  the  melting  point.  In  silicon,  this 
produces  a  f ive-fold  reduction  in  the  thermal  conductivity.  Thus 
a  nonlinear  heat  conduction  problem  must  be  solved.  We  have  also 
provided  an  analytic  solution  to  this  important  nonl inear 
problem. f.  This  problem  was  again  not  anticipated  in  the  proposal 
for  this  contract. 

Optical  Invar iants. 

An  achievement  that  was  anticipated  in  the  original  proposal 
is  a  union  of  differential  geometry  and  geometrical  optics  to 
derive  the  optical  invariants  associated  with  an  anisotropic 
crystal.  In  addition  the  surface  between  two  media  (possibly 
both  anisotropic)  has  been  treated  as  an  optical  instrument,  and 
the  full  4X4  matrix  describing  this  instrument  has  been 
obtained. 9  These  procedures  permit  the  analysis  of  brillouin 
scattering^®  and  Raman  scattering!!  experiments  needed  to  obtain 
absolute  cross-sections. 
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Ferroelectr ics. 

A  long  wave-length  Lagrangian  theory  of  the  electrodynamics 
of  pyroelectrics,  dielectrics  and  piezoelectrics  of  anv  symmetry, 
degree  of  anisotropy,  structural  complexity  w  developer  that 
describes  coupling  to  the  elastic  degrees  of  freedom  and  an  arb¬ 
itrary  number  of  internal  motions. 12  The  conservation  laws, 
boundary  conditions  and  stress  tensor  were  obtained.  Arg  sments 
based  on  momentum  conservat ion ,  stress  boundarv  condition^  and 
the  vacuum  form  of  the  Maxwell  stress  tensor  were  used  to 


establish  that  the  stress-tensor , 

not 

just  its 

d  ivergence , 

is 

meaningful,  and 

moreover  must 

be 

asvmmet r ic 

even  though 

the 

system  conserves 

angular  momentum 

13 

• 

For  a  pyroelectric,  the  boundary  conditions  at  a  moving 
surface  affect  even  the  linearized  equations  of  motion.  These 
boundary  conditions  were  obtained  bv  transforming  Maxwell's 
equations  to  the  material  frame. H 

The  linear  equations  for  elasticitv  and  piezodectr ictv  for 
a  pyroelectric  were  obtained  in  the  presence  of  a  spontaneous 
polarization  and  a  spontaneous  electric  field  Es  were 
obtained . 15 

Magnetic  Mater ials. 

The  procedures  developed  for  dielectrics  and  pvroelectr ics 
can  not  be  applied  to  magnetic  materials  without  considering  the 
following  difficulties: 
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There  are  no  obvious  variables  conjugate  to  the  three 
components  of  magnetization  (in  the  sense  that  the  mom¬ 
enta  of  the  sublattices  are  conjugate  to  the  displacem¬ 
ents  in  a  dielectric  or  ferroelectric) 


2.  The  three  components  of  magnetization  obey  the 
constraint  that  the  length  of  the  magnetization  vector 
is  fixed. 


3.  The  magnetization  is  odd  under  time  reversal  and  even 
under  inversion,  just  opposite  to  the  behavior  of  the 
polarization. 


The  constraint  (2)  has  been  handled  by  the  method  of 
Lagrange  multipl iers. The  Dirac  theory  of  singular 
Lagrangians  and  Dirac  brackets,  customarily  used  in  quantum 
field  theory,  is  then  used  to  show  that  it  is  appropriate  to 
use  the  three  components  of  magnetization  as  coordinates  and 
their  time  derivatives  as  generalized  velocities.  This 
permits  a  derivation  of  the  stress  tensor  as  well  as  the 
Lagrangian  and  Hamiltonian. 

An  extension  of  the  Green's  function  theory  developed 
for  anisotropic  d  iel ectr ics to  magnetic  materials  has  been 
developed  but  not  yet  reported.  Application  of  this  Green's 
function  and  the  fluctuation-dissipation  theorem  will  be 
used  to  calculate  Kaman  scattering  bv  magnons  or  magnaritons 
in  a  manner  analogous  to  our  previous  calculation  of  Kaman 
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scattering  by  polar itons. ** 

Detailed  applications  will  be  made  to  anti- 
ferromagnetic  as  well  as  magnetic  crystals. 

Viscous  Materials. 

The  procedures  developed  for  non-dissipative 
dielectric,  pyroelectric  and  magnetic  materials  have  been 
extended  to  viscous  materials  including  liquid  crystals. 
Our  theory  is  more  restricted  than  the  traditional, 
continuum  mechanics  approach  in  assuming  the  existence  of  a 
Lagrangian  for  the  conservative  forces  and  a  dissipation 
function  for  the  viscous  forces. 1®  This  procedure  is  more 
explicit,  and  more  deductive  than  one  based  on  the 
conservation  laws  since  the  physical  model  is  imposed  on 
construction  of  the  Lagrangian  or  dissipation  function.  Any 
approximation  made  in  the  Lagrangian,  e.g.  retention  of 
dipole  terms  but  discard  of  higher  moments,  will 
automatically  lead  to  a  consistent  set  of  truncations  in  the 
conservation  laws. 

Conservation  of  linear  and  angular  momentum  is  assured 
by  requiring  the  stored  enerqv  U  and  the  dissipation 
function  D  to  be  invariant  under  uniform  displacements  and 
rigid  body  rotations.  From  these  assumptions,  the  stress 
tensor  and  the  conservation  laws  of  momentum,  energy,  and 
angular  momentum  have  been  derived  as  well  as  the  balance  of 
entropy. 
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Initially,  the  temperature  gradient  was  used  as  an 
independent  variable  and  the  heat  flux  as  the  corresponding 
force.  The  conservation  laws  are  maintained  if  the  role  of 
these  two  variables  is  reversed. ^  However,  this  second 
choice  is  found  to  be  necessary  to  maintain  consistency  with 
the  Onsager  reciprocal  re  1  ..t  ions  .  20 

Chirp  Transform  Algorithm. 

We  have  continued  studies  of  the  electro-  magnetic 
field  distribution  inside  a  high  power  unstable  laser 
resonator.  We  have  been  mostly  concerned  of  late  with  the 
development  of  computationally  more  efficient  methods  for 
treating  edge  diffraction  at  the  output  mirror.  This  edge 
diffraction  necessitates  the  use  of  the  continuous  Fourier 
transform.  We  are  currently  developing  a  computationally 
efficient  and  flexible  continuous  Fourier  transform  based  on 
the  chirp-z  algorithm. 

Interactions  Fo_r  Multiphoton  Exc  i  ta  t  i  on . 

In  connection  with  the  coherent  multiphoton  excitation 
of  SF’g  we  have  completely  enumerated  the  twenty-two 
octahedral  invariants  involving  products  of  three  normal 
coordinates  and  also  the  ninetv-two  invariants  involving 
products  of  four  normal  coordinates  of  the  molecule.  From 
these  invariants  one  mav  construct  the  most  general  cubic 
and  quartic  cont r ibut ions  to  the  anharmonic  part  of  the  pot- 
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ential  of  the  SB  ^  molecule.  We  have  derived  a  master 
equation  describing  the  interaction  between  the  optically 
active  (/'■j  modes  in  the  presence  of  a  near-resonant  classical 
external  field.  The  anharmonic  potential  energy  (which  is 
responsible  for  the  so-called  " intra-molecul ar  relaxation") 
would  be  used  in  numerical  studies  of  the  master  equation. 
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